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Abstract 

^ . As a direct source of information on chiral symmetry breaking within 

QCD, the sigma commutator is of considerable importance. With recent 
advances in the calculation of hadron masses within full QCD it is of 
interest to see whether the sigma commutator can be calculated directly 
' from the dependence of the nucleon mass on the input quark mass. We 

, show that provided the correct chiral behaviour of QCD is respected in 

' the extrapolation to realistic quark masses one can indeed obtain a fairly 

reliable determination of the sigma commutator using present lattice data. 
, Within two-flavour, dynamical-fermion QCD the value obtained lies in the 

range 45 to 55 MeV. 
P 

' In the quest to understand hadron structure within QCD, small violations 

(~| . of fundamental symmetries play a vital role. The sigma commutator, aN- 



^N = l{N\[Q,5dQ^5,n]]\N) (1) 

■ (with Qi5 the two-flavour {i=l, 2, 3) axial charge) is an extremely important 

example. Because Qis commutes with the QCD Hamiltonian in the chiral SU(2) 
Hmit, the effect of the double commutator is to pick out the light quark mass 
term from Ti.: 

crjv = {N\ {rriuuu + m,ddd) \N) (2) 

Neglecting the very small effect of the u — d mass difference we can write Eq. (||) 
in the form 
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CTN = {N\fa{uu + dd)\N) (3) 
dm 



"^^^ (4) 
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with m = (ttIu + mdj/'i. Equation (Q) follows from the Feynman-Hellman the- 
orem 1^. 

While there is no direct experimental measurement of ctat, the value inferred 
from world data has been 45 ± 8 MeV |^ for some time. Recently there has 
been considerable interest in this value because of progress in the determination 
of the pion-nucleon scattering lengths ^ and new phase shift analyses ||^, ^ . 
For an excellent summary of the sources of the proposed variations and the 
disagreements between various investigators we refer to the excellent review of 
Kneckt For our purposes the experimental value is of limited interest as 
the full lattice QCD calculations upon which our work is based involve only two 
active flavours. Nevertheless, as a guide, the current work suggests that the 
best value oi gn may be between 8 and 26 MeV larger than the value quoted 
above |^. 

Numerous calculations of (Tjv have been made within QCD motivated models 
[0 and there has been considerable work within the framework of chiral pertur- 
bation theory |Q. However, direct calculations of within QCD itself have 
proven to be difficult. Early attempts |[l2) to extract (Tat from the quark mass 
dependence of the nucleon mass in quenched QCD (using Eq.(^ ) produced 
values in the range 15 to 25 MeV. Attention subsequently turned to determin- 
ing ctat by calculating the scalar matrix element of the nucleon {N\uu + dd\N). 
There it was discovered that the sea quark loops make a dominant contribution 
to (Tat These works, based on quenched QCD simulation, found values 

in the 40 to 60 MeV range, which are more compatible with the experimental 
values quoted earlier. 

On the other hand, the most recent estimate of ctjv, and the only one based 
on a two-flavour, dynamical-fermion lattice QCD calculation, comes from the 
SESAM collaboration. They obtain a value of 18 ± 5 MeV @, through a 
direct calculation of the scalar matrix element {N\uu + dd\N). The discrepancy 
from the quenched results of Refs. [l5| is not so much an unquenching 
effect in the scalar matrix element but rather a signiflcant suppression of the 
quark mass in going from quenched to full QCD. The difficulty in all approaches 
which evaluate {N\uu + dd\N) is that neither it nor m is renormalization group 
invariant. One must reconstruct the scale invariant result from the product of 
the scale dependent matrix element and the scale dependent quark masses. The 
latter are extremely difficult to determine precisely and are the chief source of 
uncertainty in this approach. 

An additional difficulty in extracting aN from lattice studies is the need to 
extrapolate from quite large pion masses, typically above 500 or 600 MeV. An 
important innovation adopted by Dong et al., but not by the SESAM collabo- 
ration, was to extrapolate the computed values of {N\uu + dd\N) using a form 
motivated by chiral symmetry, namely a + bfh^ . On the other hand, the value 
of b used was not constrained by chiral symmetry and higher order terms of the 
chiral expansion were not considered. Furthermore, since the work was based 
on a quenched calculation, the chiral behaviour implicit in the lattice results 
involves incorrect chiral coefficients | [T^ . 
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Our work is motivated by recent, dramatic improvements in computing 
power which, together with the development of improved actions mean 
that we now have accurate calculations of the mass of the nucleon within full 
QCD (for two flavours) as a function of m down to to^ 500 MeV. (Since 
is proportional to to over the range studied we choose to display all results as 
a function of mj.) In addition, CP-PACS has recently published a result at 
TOtt ~ 300 MeV , albeit with somewhat large errors. Provided that one has 
control over the extrapolation of this lattice data to the physical pion mass, 
TO^ = M = 140 MeV, one can calculate (Jn by evaluating Eq. (Q) at the physical 
pion mass. Note that this approach has the important advantage over the cal- 
culation of the scalar density that one only needs to work with renormalization 
group invariant quantities. We therefore turn to a consideration of the method 
of extrapolation. 

The lattice data for the nucleon mass calculated by UKQCD ^ and CP- 
PACS dl is shown in Fig. |l|. Both groups cite a 10% uncertainty in setting the 
lattice scale, so we have scaled the former down and the latter up by 5% so that 
the data sets are consistent. Over almost the entire range of toJ, the data shows 
a dependence on quark mass that is essentially linear. However, the preliminary 
point at 0.1 GeV^ suggests some curvature in the low mass region. This is 

indeed expected on the basis of chiral symmetry with the leading non-analytic 
(LNA) correction (in terms of to) being proportional to m% (m'^/^): 

<5Af-* = 7-*m3 , = . (5) 

These observations led the CP-PACS group to extrapolate their data with the 
simple, phenomenological form: 

Mn + (3ml + ^ml . (6) 

The corresponding fit to the combined data set, using Eq. (||), is shown as the 

short-dashed curve in Fig. |^ and the parameters (a, /3, 7) = (0.912, 1.69, —0.761) 

(the units are appropriate powers of GeV) . This yields a value for the sigma com- 
(p) 

mutator, apf = 29.7 MeV, where the superscript stands for "phenomenological". 

The difficulty with this purely phenomenological analysis was discussed in 
Ref. That is, the value of 7 = —0.761 is almost an order of magnitude 

smaller than the model independent LNA term, 7"""* = —5.60 GeV~^. Clearly 
this presents some concern when evaluating a^, because of the derivative re- 
quired. An alternative approach to this problem was recently suggested by 
Leinweber et al. |[l8). They realised that the pion loop diagrams. Fig. ||(a) and 
||(b) not only yield the most important non-analytic structure in the expression 
for the nucleon mass, but amongst all the possible meson baryon states which 
contribute to the nucleon mass within QCD, they alone give rise to a significant 
variation of the nucleon mass as m^r ^ 0. In Ref. ||l8j it was suggested that one 
should extrapolate Mn as a function of quark mass using: 

Mn = a + /3to^ + ctat at (to^ , A) + ctat a (tott , A) , (7) 
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Figure 1: Nucleon mass calculated by CP-PACS (solid points) and UKQCD 
(open points), as a function of mj, both are scaled by 5% to improve consistency. 
The soHd curve is a fit to Eq. (^ with a 1.225 GeV dipole form factor, while the 
dashed curve is the same fit using a sharp cut-off form factor (^?-function) . The 
short-dash curve is a fit to Eq. (^, and the long-dash curve is a fit to Eq. (^ 
excluding the lowest data point. The vertical line indicates the physical pion 
mass. 



where (tnn and ctjva are the self-energy contributions of Figs. y(a) and |2|(b), 
respectively, using a sharp cut-off in momentum, 9{A — k). The full analytical 
expressions for aMN and ctata are given in Ref. |Q. For our purposes it suf- 
fices that they have precisely the correct LNA and next-to-leading non-analytic 
behaviour required by chiral perturbation theory as m„ 0. In addition, ctata 
contains the correct, square root branch point (~ [mj — (Ma — M^Y]'^) at the 
A — N threshold, which is essential for extrapolations from above the A — N 
threshold. 

Fitting Eq. (|^) to the data, including the point near 0.1 GeV^, gives the dot- 
dash curve in Fig. |^ ((a,/3. A) = (1.42,0.564,0.661)). The corresponding value 
of aN is 54.6 MeV and the physical nucleon mass is 870 MeV. Omitting the 
lowest data point from the fit yields the long-dash curve in Fig. |l| ((a,/3. A) = 
(1.76, 0.386, 0.789)) with — 65.8 MeV. Clearly the curvature associated with 
the chiral corrections at low quark mass is extremely important in the evaluation 

of (TJV. 

In order to estimate the error in the extracted value of cr at we would need to 
have the full data set on a configuration by configuration basis. As this is not 
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Figure 2: One-loop pion induced self energy of the nucleon. 



Scaling 


(TJV 


CP-PACS 


UKQCD 


Dipole 


Sharp 


Cubic 


5 


-5 


47.2 ±1.8 


54.6 ±2.0 


29.7 


10 





48.1 ±1.9 


54.9 ±2.0 


28.6 





-10 


45.4 ±1.9 


54.3 ±1.9 


31.0 



Table 1: Sigma Commutator Values. The Dipole and Sharp results were 
calculated with our preferred form of a + /3m^ + o'ArAr(A, rrtTr) + o'ArA(A, m^r) 
with either a dipole form factor for the Ntt vertex or a 6'-function. The values 
of dipole parameter (A^i) were (1.225, 1.250, 1.175) GeV. The Cubic results are 
for the a + I3rn^ + 7771^ extrapolation function, with 7 unconstrained by chiral 
symmetry - as explained in the text this produces an unreliable value for cfn ■ 



available, the errors that we quote are naive estimates only. The extracted value 
of gn is very well determined by the present data, the result being 54.6 ± 2.0 
MeV. Since the process of setting the physical mass scale via the string tension 
is thought to have a systematic error of 10%, one might naively expect this to 
apply to (TAT. However, all masses in the problem including the pion (or quark) 
mass, as well as that of the nucleon, scale with the lattice parameter a. It turns 
out that when one uses Eq. (^ at the physical pion mass (which means a slightly 
different value of fha if a changes), the value of un is extremely stable. If, for 
example, one raises the CP-PACS data by 15% and the UKQCD data by 5% 
(instead of 5% and —5%, respectively) the value of (Jn shifts from 54.6 ± 2.0 to 
55.2 ±2.1 MeV. We present calculations in Table ^ that show, for a variety of 
scaHngs of the lattice data, how stable our results are. 

The remaining issue, for the present data, is the model dependence associated 
with the choice of a sharp cut-off in the pionic self-energies. Our investigations 
in Ref. ||l^ showed that Eq. (^ could reproduce the dependence of Mn on 
within the cloudy bag model, and that it could also describe the dependence of 
pion self-energy terms calculated with dipole form factors. Thus we believe that 
any model satisfying the essential chiral constraints and fitting the lattice data 
should give essentially the same answer. We checked this by numerically fitting 
the lattice data (solid curve) with the form of Eq. ^ but with onn and (Tata 
calculated with dipole form factors of mass A^) at all pion-baryon vertices. Since 
the preferred phenomenological form of the Ni: form factor is a dipole, we regard 
the dipole result shown in the first line of Table [| as our best estimate, namely 
aN = 47.2 ± 1.8 MeV with fit parameters (a,/3,A£,) = (2.02,0.398,1.225). A 
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remaining source of error is that, although the lattice results were calculated 
with an improved action, there still is an error associated with the extrapolation 
to the infinite volume, continuum limit. The importance of the inclusion of the 
correct chiral behaviour is clearly seen by the fact that it increases the value 
of the sigma commutator from the 30 MeV of the unconstrained cubic fit to 
around 50 MeV. 

Clearly an enormous amount of work remains to be done before we will fully 
understand the structure of the nucleon within QCD. It is vital that the rapid 
progress on improved actions and faster computers continue and that we have 
three fiavour calculations within full QCD at masses as close as possible to the 
physical quark masses. Nevertheless, it is a remarkable result that the present 
lattice data for dynamical-fermion, two-flavour QCD, yields such a stable and 
accurate answer for the sigma commutator, an answer which is already within 
the range of the experimental values. The implications of this result for models 
of hadron structure need to be explored urgently. 

One of us (SVW) would like to acknowledge helpful discussions with Tom 
Cohen at an early stage of this work. We would also like acknowledge helpful 
comments from Chris Allton, Craig Roberts and Robert Perry. This work was 
supported in part by the Australian Research Council. 
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